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1. The exponential convexity result 

Definition. A function / on K, / : R —» [0, oo), is said to be exponentially 
convex if 

1. For every positive integer N, for every choice of real numbers t±,t 2 , ■ ■ ■, 
t]y, and complex numbers £i, £ 2 , • ■ • the inequality holds 

N 

2 f(t r + t s )trZ> 0; (1.1) 

r,s= 1 

2. The function f is continuous on R. 

The class of exponentially convex functions was introduced by S.N.Bernstein, 
[B], see §15 there. Russian translation of the paper 0 can be found in EH 
pp.370-425]. 

From (11.111 it follows that the inequality /(ti + ^ 2 ) ^ 'y//(2ti)/(2t2) 
holds for every t\ e R, £2 6 K- Thus the alternative takes place: 

If f is an exponentially convex function then either f(t) = 0, or f{t) > 0 for 
every f e R. 
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Properties of the class of exponentially convex functions. 

PI. If /(f if an exponentially convex function and c ^ 0 is a nonnegative 
constant, then the function cf(t) is exponentially convex. 

P2. If fi(t) and / 2 (f) are exponentially convex functions, then their sum 
fi(t ) + / 2 (f) is exponentially convex. 

P 3. If /i(t) and / 2 (f) are exponentially convex functions, then their product 
fi(t) ' f‘ 2 (t) is exponentially convex. 

P4. Let {fn{t)}i^ n <oo be a sequence of exponentially convex functions. We 
assume that for each leR there exists the limit f(t) = lim.n-.oo / n (t), 
and that f(t) < oo Vt G R. Then the limiting function f(t) is exponen¬ 
tially convex. 


From the functional equation for the exponential function it follows that 
for each real number A, for every choice of real numbers t±,t 2 , ■ ■ ■, ijv and 
complex numbers £i, £ 2 , • ■ • , £jvj the equality holds 

N N 2 


2 


g\(t r +t s ) 





( 1 . 2 ) 


r,s= 1 


P=1 


The relation (11.21) can be formulated as 


Lemma 1.1. For each real A, the function e xt of the variable t is exponentially 
convex. 


For z G C, the function coshz, which is called the hyperbolic cosine of 
z, is defined as 

cosh z = i (e z + e~ z ). (1.3) 

From Lemma 1 1.1 1 and property P2 we obtain 

Lemma 1.2. For each real fi, the function cosh (/it) of the variable t is expo¬ 
nentially convex. 

The following result is well known. 

Theorem (The representation theorem). 

1. Let cr(d\) be a nonnegative measure on the real axis, and let the function 
f(t) be a two-sided Laplace transform of the measure dcr(X): 

f{t)= j e xt da{ A) (1.4) 

AeR 

for any f el. Then the function f is exponentially convex. 

2. Let f{t) be an exponentially convex function. Then this function f can be 
represented on R as a two-sided Laplace transform of a nonnegative 

measure d<r{X). (In particular, the integral in the right hand side of (11.41) 
is finite for any t G R.) The representing measure dcr(X) is unique. 

The assertion 1 of the representation theorem is an evident consequence 
of Lemma o of the properties P1,P2, P4, and of the definition of the 
integration. 




The function cosh (fa t 2 + b'j is exponentially convex. 
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The proof of the assertion 2 can be found in m , Theorem 5.5.4, and in 
[Wi] . Theorem 21. 

Of course, Lemma fL2l is a special case of the representation theorem 
which corresponds to the representing measure er(dA) = 1/2(<5(A — /x) + <5(A + 
/ i )) dX, where 6( A + p.) are Dirak’s ^-functions supported at the points +p. 


Lemma 1.3. The expression 

y(f, a , b) = cosh (V at 2 + b ). (1.5) 

is well defined for every complex numbers t, a, b. The function ip(t , a, b) is an 
entire function of complex variables ( t , a, b) e C 3 . For each fixed a > 0 and b, 
the function tp(t,a,b), considered as a function oft, is an entire function of 
exponential type yet. 

Proof. The function ip(t, a , b) is a superposition of the entire function cosh 
of variable ( and the quadratic polynomial </( t , a, b) = at 2 + b. The assertion 
concerning the growth of this function is evident. □ 

In the paper [BM VI a conjecture was formulated which now is commonly 
known as the BMV conjecture: 

The BMV Conjecture. Let A and B be Hermitian matrices of size n x n. 
Then the function 

fA,B(t) = trace {exp[tH + B]} (1.6) 

of the variable t is representable as a bilateral Laplace transform of a non¬ 
negative measure daA,B{ A) compactly supported on the real axis: 

fA,B(t)= j exp(iA) da A ,B(X), Vte (-00,00). (1.7) 

Ae(—00,00) 

In general case, if the matrices A and B do not commute, the BMV con¬ 
jecture remained an open question for longer than 35 years. In 2011, Herbert 
Stahl, [St], gave an affirmative answer to the BMV conjecture. 

Theorem (H. Stahl) Let A and B be nx n Hermitian matrices. Then the func¬ 
tion f a b ( t ) defined by (11.61) is representable as the bilateral Laplace transform 
G2D of a non-negative measure dcrA,B(A) supported on the closed interval 

[A m i n , A max ] - 

The proof of Herbert Stahl is based on ingenious considerations related 
to Riemann surfaces of algebraic functions. 

In the case 2x2 matrices A and B , the BMV conjecture is equivalent 
to the exponential convexity with respect to t for each a > 0 and 6 > 0 of 
the function ip(t, a , b ) = cosh (V at 2 + 6) which was introduced in (L5]) . 

Theorem 1. For each fixed a > 0 and b ^ 0, the function y(t, a, b) defined by 
C3D is an exponentially convex function of variable t. 

In what follow we present three different proofs of Theorem [0 The first 
and the second proofs are based on the representation theorem. We prove 
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that the function d( A, a, b) which is defined by (11.81) below takes positive 
values for A e (—y/a, V®)- I n the first proof we calculate the function d( A, a, b) 
explicitly expressing this function in terms of the modified Bessel function I\. 
In the second proof, we prove the positivity of the function d{ A, a, b) using 
the reasoning by Herbert Stahl in [St]. (We use a very simple special case 
of this reasoning.) The third proof is based on the Taylor expansion of the 
function ip(t,a,b ), (11.51) . with respect to parameter b. This proof does not 
use any integration in the complex plane. It based only on Lemma 11.21 and 
on the properties PI - P4 of the class of exponentially convex functions. 
As a by-product of this proof we obtain that all coefficients of this Taylor 
expansion are exponentially convex functions. However we can not conclude 
directly from this proof that the restriction of the representing measure on 
the open interval yfa) is an absolutely continuous measure. 

Lemma 1.4. For each fixed a > 0, b ^ 0, the function ip(t, a, b ) defined by 
(O) is representable in the form 

•Jo, 

<p(t, a, b) = cosh yfa, t + / d(X,a,b)e xt dX, Vie C, (1.8) 

— -Ja 

where the function d{ A, a, b) possesses the properties 


1 . 

2 . 

3. 


rJa 


I —•Ja 


|d(A, a, b)\ 2 dX < go; 


(1.9) 


The function d( A, a, b ) is continuous with respect to X on the closed 
interval [—y/a,y/a\, takes real values there, and is even. 

The values of the function d( A, a, b ) at the end points + y/a of the interval 
[--v/a, Va] are: 


d(+y/a, a , b) = 



( 1 . 10 ) 


Proof. We introduce the function 

d(t, a, b) = cosh (V at 2 + b) — cosh \fat (1.11) 

of variables t, a, b. Considered as a function of t for fixed positive a and b , 
d(t, a, b) is entire function of exponential type \fa. On the imaginary axis d 
takes the form 

d(ir, a, b) = cos V ar 2 — b — cos Var, re R. (1.12) 

From (] 1. 121) it follows that the function d is a bounded and decaying on the 
imaginary axis: |d(zr, a, b) ^ 1 + coshfr, r e R, d(ir, a, b) = 0(|r| _1 ) as r —> 
±oo. By the Wiener-Paley theorem, the function d(ir,a,b) is representable 

in the form _ 

•Ja 

d(ir,a,b)= j d(X,a,b)e lXr dX, rel, (1.13) 

—yfa 







The function cosh ( yfat 2 + b'j is exponentially convex. 
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where the function d( A, a, b) satisfies the condition (11.91) . The equality (11.131) 
serves as a definition of the function d(X,a,b). So, this function is defined 
only for a > 0, b > 0, —<Ja < A =$ y/a. 

Since the function d(ir, a , b ) is even with respect to r and real valued, its 
inverse Fourier transform d(A, a, b ) is even with respect to A and real valued. 

From (11.121) we obtain that 


d{ir , a, b) 


b sin y/a t 
2 y/ar 


0(t 2 ) as t 


+ 00 . 


(1.14) 


Hence the function in the left hand side of (11.1411 is a Fourier transform of 
some function ?'(A) which is square summable and continuous at every A e R. 
We remark that 


b sinY^r 
2 y/ar 


\fa 



—y/a 


re R, 


where c = -/t 
4 y/CL 


is a constant function. Hence 


r-m = J d ( A ’ a ’ fo ) 

jo 

Since r(A) = 0 for |A| > y/a, also r(± v / a) 


for |A| < y/a, 
for |A| > y/a. 

0. Thus, (11.101) holds. 


□ 


2. Representation the function d( A, a, b ) 
by a contour integral. 

Let S' be a segment of the imaginary axis: 

S={c = £ + «?: e = 0,-^1 < < V? }• (2 - 1} 

The function yjaC/ + b is a single value function of ( in the complex plane 
slitted along the vertical segment S. We choose the branch of this function 
which takes positive values for large real (. 

Lemma 2.1. The function d(X,a,b), which was defined by (11.131) . admits the 
integral representation 

d(X,a,b) = -— (^e - '^ a< ’ 2+fc e -Al ’ —y[a<X<y[a 1 (2.2) 

47T'j J 
r 

where T is an arbitrary counterclockwise oriented closed Jordan curve which 
contains the slit S inside. 

Proof. According the inversion formula for the Fourier transform, 

00 

d(X,a,b) = ^ J (p(ii), a, b)e~ lXri dip 
— 00 


(2.3) 
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We interpret the integral in the right hand side of (12.31) as the integral along 
the the vertical straight line {C : Re£ = 0}: 

+iR 

d(\, a, b) = —[ ip(C,a,b)e~ x<: d( = lira f ip{(,a,b)e~ xc d(. 

27 n J R^+co 2m J 


ReC=0 


— iR 


(2.4) 

Since the function ip((, a , b) is bounded in each vertical strip {C : a < Re C ^ 
/3} and tends to zero as ImC —> ±oo within this strip, the value of the integral 
in (12.41) does not change if we integrate along any vertical line {C : Re C = 7 }, 
where 7 is an arbitrary real number: 

d(\,a,b) = -—: [ <£>(£, a, b)e~ x ^ dC,, — Va^A^v/a. (2.5) 

2m J 


ReC=7 


Choosing 7 < 0 , we split the integral in (12.51) into the sum 

d(A, a, 6 ) = ^7 J <p+((,a,b)e~ Xi d( + J <p-{C,a,b)t 


-AC 


dC, 


Re C=7 


Re C = 7 


( 2 . 6 ) 


where 


<p-(C,a,b) = 


_ 1 ( p -\/a-i 2 +b _ „-VaC 


)• 


(2.7) 

The function ip + ((,a, b) is holomorphic in the halfplane {C : Re^ < 7 } 
and admits the estimate 

l¥>+(C,M)l < c( 7 )(l + |C|)“ 1 e V “ ReC , VC : ReC < 7 
there, where c( 7 ) < 00 is a constant. Therefore 

7 ^- J ^+(C,a,b)e~ x( d( = 0, A < ■s/d. (2.8) 

Re C = 7 

The function </?_ (C, a , b) is holomorphic in the slitted half plane {C : Re C S 5 7 , 
C £ S} and admits the estimate 

|^_(C,a,&)| c( 7 )(l + |C|) _ 1 e _V “ ReC , VC : ReC > 7, C $ $ 

there. Therefore 

J <P-((,a,b)e- x<: d( =(p-((,a,b)e~ x<: e- xc dC, A >-y/a, 

Re C = 7 r 

(2- 9 ) 

where T is an arbitrary closed Jourdan curve which is oriented clockwise 
and contains the slit S in its interior. Since the function is entire, 

$ e-V^C = 0 . So 
r 

a, b)e~ x< ^ dQ = -^- ^ e -V^> e -*C dC, (2.10) 





The function cosh (V at 2 + 6) is exponentially convex. 
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where the integral in the right hand side is taken over the curve T which is 
oriented counterclockwise. Comparing (12.61) , (12.81) , (12.91) , and (12.101) , we obtain 

(EUD- □ 


3. Explicite calculation of the function d( A, a, b). 

Lemma 3.1. The function d(X,a,b) which was defined by (11.131) admits the 
integral representation 

d(X, a, b) = — I ^ J sinh \/b(l — r 2 ) • cos dr, —\[a ^ A ^ yfa . 

o 

(3.1) 

Proof. We derive Lemma 13.11 from Lemma 12.11 showing that 


- ^-<^e-V“C 2 + b e - A <: d( = lyJlL J sinhVKl-^ 2 ) • cos A yf^rdr. 
r o 

(3.2) 

The function e~^ a ^ 2+b e~ x< ^ is holomorphic in the domain C\S and contin¬ 
uous up to boundary S = d(C\S) of this domain. Therefore the integral of 
this function over T does not change if we shrink the original contour T to 
the boundary S: 

_L I e -^ a C 2 + b e -K d(= (^e-V“C 2 + b e - A C d( (3.3) 

47T i J 47rj Jr 

r s 

To one "geometric" point ip e S there corresponds two topologically different 
"boundary" points +0 + ip and —0 + ip lying on the right edge S + and the 
left edge S~ of the slit S respectively. The chosen branch of the function 
yjaf 2 + b takes the following values on the boundary of the domain C \S: 

yj a(+0 + ip) 2 + b = — \J a(—0 + ip) 2 + b = y/b — ap 2 , ip e S. (3.4) 


If the point £ = ±0 + ip runs over S = d ( C\S ) counterclockwise, the p 


increases from 
f =e S~. Therefore 


to y A if £ =e S + and p decreases from 


to 


if 


yfbfa, 


j' e~^+~ b e" AC d( = + i J e -^ b - ar,2 e- lXv dp, 


s+ 


-y/b/a 

y/bfa. 


j- e -V“C 2 + b e " AC d( = -i J e + ^ b ~ ari2 e- iXv dp 

S~ -y/b/a 
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Thus 


y/bfa 

— <k e -V“C 2 + b e - AC d( = — [ (e-V 6 -^ 2 - e s/ b ~ a v 2 \ e -i*v dr] = 

TT'i J 47 T J 

-y/b/a 


y/bja yfbja 

— / sinh yfb — arj 2 e~ lXri dr] = -/ sinh yfb — arf cos A rj dr] = 

27 t J irj 

-yfiR 0 

1 

J sinh^/fe(l — d 2 ) ■ c ° s (Vi^) ^ ^ 3 - 5 ) 


Comparing (12.21) with (13.51) . we obtain (EH). 

Lemma 3.2. Let a > 0 and b > 0 be fixed positive numbers. Then 


1. The function d(X,a,b) which was defined by (11.131) can be expressed 
explicitly in terms of the modified Bessel function I\: 


g?(A, a, b) = 


yfb 


( a—\ 2 )b 


-yfa < A < yfa. (3.6) 


2 y/a - A 2 

2. The function d{ A, a, b ) is representable as the sum of the series 
b fX 1 /(a — A 2 )6 


d(A ’ a,&) 4VS fe ? 0 fc!(fc + l)! 


4a 


> 0, — Va < A < yfa, b ^ 0. (3.7) 


Remark 3.1. TTie expression in the right hand sides of (EH is an entire 
function of three variables (A, Va -1 , b) e C 3 . However the equalities (13.11) . 
m , eh hold only for a > 0, b > 0, — yfa ^ A < yfa. (We recall that the 
function d(X,a,b) was defined by (11.131) only for a > 0, b > 0, —yfa < A 
yfa.) 


Proof of Lemma 13.21 We start from the formula (EH- Using the Taylor ex¬ 
pansion of the hyperbolic sinh function, we obtain 


d{ A, a, b) = 


y- 1 

£ n (2r+l)! 


*/ (1 


— r 


2 A 7 -H 


cos 


(Vm 


dr (3.8) 


7r 
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The integral in the right hand side of (13.81) can be expressed in terms of the 
Bessel function J r +i, see [Ml 9.1.20]: 


J (l-r 2 ) r+ 5cos(Ay^r) dr 


= 7r 1 / 2 2 r T(r+ 3/2)a r ^ i 6- £ f L ) A- (r+1) J r+1 (\^j. 
Substituting the last equality into (13.81) . we obtain the equality 
_I W iXr + 3/2) r r±l / 

= 7T 2 \ -1- L -19'0 2/1 2 X t 


d(\a,b) = n * £ -r. + - y a* 

r=0 A ' 




Taking into account the duplication formula for the Gamma-function, 

6.1.18]: 

T(r+|) = 7j _i 2 _ (2r+1) _ 1 


T(2r + 2) T(r + 1) ’ 

we transform the last equality to the form 

00 1 /— 

d(X,a,b) = £ -2-^aib^ L \~^J r+l (\^^. (3.9) 


r=0 


Now we would like to reduce the equality (|3.9D to the form which occurs in 
the so called Multiplication Theoren 0, see |AS[ 9.1.74]: 


aa)-#| 


r=0 


■ ifVf ) ^ + i(A 


Let us introduce p : p 2 — 1 = — —, i.e. 


A 2 ’ 


/i = i- 


.Va — A 2 


A 


(3.10) 


(3.11) 


Then the equality (13.101) takes the form 

S{Ka - b) - w=b? ■ '‘fU*'* 2 - 'I’tiVi)’*«(Vi)- (3-i2) 

According to the Multiplication Theorem, 

1 
r 


^h ( ^ 2 -!) r 0 Vl) r J v + i( aV!) = (*V^a^v^) > VA ’ a ’ 


r=0 


(3.13) 

Taking into account that J\{iz) = Hi{z) : we reduce the equality (13.121) to 
the form (13.61) . 


1 Proof of the Multiplication Theorem can be found in no Chapter V, sec. 5.22], see formula 

(15) on page 142 of the English edition or on the page 156 of the Russian translation. See 
also IScl Chapter IV, sec.21]. 

















10 


V. Katsnelson 


Using the Taylor expansion of the modified Bessel function I\ , (A~S , 
9.6.10], we represent the function d( A, a, b) as the sum of the series (Id.71) . □ 


4. The first proof of Theorem [B 

From the equality lETTll is evident that 

d( A, a, b) > 0 for A e [—Va, V^] (4-1) 

Theorem [l] follows from (14.11) and (11.81) . □ 


Theorem 2. For each a > 0, the function ip(t,a,b) which was introduced in 
admits the Taylor expansion with respect to b: 


o° 1 

<£>(*,a, b) = s — ipk(t,a)b k , VteR, 0 =% 6 < go. 
k =o 


(4.2a) 


For each k ^ 0, the function ipk(t,a), which is the k-th coefficient of the 
Taylor expansion (I4.2all . is exponentially convex: 

ipo(t,a) =coshy / at, (4.2b) 


Tk(t,a) = 


(fc — l)!4 fc a‘ 


\J~a- 

— f (a — A 2 ) 
■ 2 J 


k —1 


dA, 


—yfa 


k= 1,2,3,.... (4.2c) 


Proof. The expansion (13.71) can be presented as a Taylor expansion with 
respect to b: 


d(X,a,b) = £ ld fc (A,a)& fc , 


fe=i 


(4.3a) 


where 

dfc(A.a) = ( °~*2 fc \ , fc = 1,2,3,.... (4.3b) 

(A: — l)!4 /c a 2 

Substituting the expansion (|4.3[) into the integrand in (|1.8|) . we obtain the 
expansion (ixa . It is evident that c4(A,a) > 0 for —y/a < A < -^/a. The 
exponential convexity of the function ip k (i, a) follows from the representation 

m . □ 


Remark 4.1. The function ip k {t,a), (|4.2c|l . can 6e expressed in terms of the 
modified Bessel function I k _ 1 : 

a) = 7r 5 2 _ ^ fe+3 ^a _ ^ +3 t _ ^ fc_5 ^ J fc _i (Vaf), a > 0, t e R, k = 1,2,3,.... 

(4.4) 


See ESI 9.6.18], 










The function cosh (\/at 2 + b) is exponentially convex. 11 

Remark 4.2. The formula (13.71) appeared in [StI subsection 7.3], see formulas 
(7.22) and (7.23) there. In [St], the derivation of the expansion (13.71) was done 
by a direct calculation, without any reference to the multiplication theorem for 
Bessel function. It should be mention that the series in the right hand side of 
433 appeared in jMKi section 2] as a perturbation series related to the BMV 
conjection for 2x2 matrices. 


5. The second proof of Theorem [B 

The starting point of the first as well as of the second is the representation 
of the value d{ A, a, b) by the contour integral (12.21) . See Lemma Im 

In the first proof, we shrank the contour of integration over the slit S, 
so the contour of integration was the same for every A e [— yfa, v / a]. 

In contrast to this, in the second proof we choose the contour T in such 
a way that the exponent — y/ aCf 2 + b — A£ of the integrand in 

(]2.2I) takes real values on T. (So the contour T depends on A!). We denote 
this contour by Ta 

The function d(A, a , b ) is even with respect to A. Therefore to prove the 
exponential convexity of the function ip(t, a , b), it is enough to prove that the 
value g?(A, a, b ) is positive for each 

Ae(-Va,0). (5.1) 

We choose an arbitrary A satisfying the condition (15.11) and fix this choice in 
the course of the proof. 

Let us introduce the functions 

u(C) = Re (\/af 2 + b + AC), z e C\S , (5.2a) 

u(C) = Im(VaC 2 + b + AC), z e C\S, (5.2b) 

where S is the vertical slit (ED and the branch of the function y/ af 2 + b in 
C\S is chosen which takes positive values for large real C- 

Lemma 5.1. Let us assume that a > 0, b > 0 and A satisfies the condition 
(ED- Then there exist e > 0 R < oo, e = e(a, b , A), R = R(a , b , A), such that 

u(C)/ImC>0, VC e C : |C| > R, ImC ¥= 0, (5.3a) 

u(C)/InrC < 0, VC e C\S : |C| < e, ImC ^ 0. (5.3b) 

Proof. From the identity 


+ b 

we derive that 



- VaC = 


yjaf 2 + b + VaC 


Im a( 2 + b — y/a Im C = — (Im y/af 2 + b + yfa Im C) p(C) 
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where p(Q = b\yf aCf 2 + b + y/at\ 2 . Thus 


and 


Im yjaC , 2 + b = 1 ■ yfa Im C 

! + P(C) 

u(£) = (Va 7——777 + a') • ImC- 


V i + P(C)'-7— (5 - 4) 

It is clear that@ p(£) —» 0 as |£| —» 00, p(() —» 1 as |£| —» 0, ( f S. Since 
yfa + A > 0, the inequality (1 5. dal) holds if if |£| is large enough. Since A < 0, 
the inequality (I5.3bl) holds if |£| is small enough. □ 

Let N\ be the set 


N x = {( e C\S : v(() = 0}, 

Lemma 5.2. 

1. The set N\ is the union of the real axis and an ellipse P A : 

Nx = (K\0) + r A . 

where the ellipse P A is described by the equation: 


^ + W = 1, (C = <£ + «?), 


with 


A = \ l — ■ J 4 L( 1 -- 


B = a - • 1 - 


a yfa 

2. The slit S is contained in the interior of the ellipse T A . 

Proof. 


(5.5) 

(5.6) 

(5.7) 

(5.8) 


1. Let £ = £ + ip, yjaC, 2 + b = p + iq, where £, r],p, q are real numbers. The 
equality 

±V a £ 2 + b = p + iq 
is equivalent to the system of equalities 

\a(Z 2 -V 2 ) + b = p 2 ~ q 2 , 

[ a£r) = pq. 

Here p = pf, rj), q = q(£, ry). Clearly rr(£, ry) = qf, ry) + Ary. 

Let £ e N\. This means that r;(£,ry) = 0, i.e. 

q = —Ary (5.10a) 

Substituting this equality into the second equality of the system m , we 
obtain the equality a£ry = —A pr). Assuming that ry ^ 0, that is £ f K, we can 
cancel by ry and obtain 

at, 


P = -~ 


(5.10b) 


2 Here the choice of the branch of the function yj a£ 2 + b is important. 
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Substituting the equalities (15.101) into the first equality of the system (15.91) . 
we obtain that the equality (EH holds for ( = £ + ip. Thus we proved that 

(N X \R) c T A . (5.11) 

Let 

H+= {C : ImC > 0}, = {C : ImC < 0} (5.12) 

be the upper and the lower half-plane respectively. 

According to Lemma f5. 11 there exist points £ 6 H + \5' where v(Q > 0 
and points ( G H + \S' where v(() < 0. This means that the set N\, (15.51) . 
separates the domain EI + \S'. In other words, the open set (H + \iS , )\./V a is dis¬ 
connected. Since v(f) = —v((), the set N\ is symmetric with respect to the 
real axis. The set T A \1V A also is symmetric with respect to the real axis. Since 
(15.111) . the set N\ can not separate the domain (H + \S I ) if T\\N\ ^ 0. 

2. In view of (EH), the inequality 0 < A < B hold. So A is the minor semi¬ 
axis of the ellips T A and B is its major semiaxis. Moreover, the inequality 

< B holds. This means that the slit S is contained inside the ellipse T A . □ 

Lemma 5.3. 

1. The functions u(f) and v(f) are conjugate harmonic function of £ in 
the domain £ G C\S. 

2. The only critical points of the the functions u and v in the domain 
£ G C\S are the points 

a ) 

(5-13) 

that is the points where the ellipse T A and the real axis R intersect. 

3. If C G H + lies outside the contour T A) then v(f) >0. If £ G H + \S' lies 
inside the contour T A) then v(Q < 0. 

Proof. The functions u and v are the real and the imaginary parts of the 
holomorphic function -\Jaf 2 + b + A£. From the Cauchy-Riemann equation it 
follows that the functions u and v have the same critical points. Moreover 
the point ( is critical for v if and only if (j is a root of the derivative af(aC, 2 + 
b)~z + A of the function af 2 + b + AC- An explicit calculation shows that 
this derivative has only two roots C+(^) an d C-(^), (15.131) . 

Let Ext(T A ) and Int(T A ) be the exterior and the exterior of the contour 
T A respectively. Each of the sets Ef xt and 

£f xt = Ext(T A ) n H+, E{ nt = Int(T A ) n (H^S) (5.14) 

is a connected open set. According to (EH and (EH), the continuous real 
valued function v does not vanish on any of these two sets. Hence the values 
v(() have the same sign, say s Ext , at all points £ of the set and the 

same same sign, say s 7nt , at all points f of the set £ A nt - Now the assertion 3 


C + (A) = 


|A| 


2 \ -2 


1 - 


and C-(A) = — 


b |A| 
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of Lemma 15.31 is a consequence of Lemma 15.11 □ 

Completion of the proof of Theorem |T| Let us chose the ellipse Ta as the 
contour of integration L in the integral in the right hand side of (12.211 . Since 
the imaginary part v(Q of the exponent of the integrand vanishes on Ta, the 
integral representation (12.21) takes the form 

d(X, a, b) = — -j— (jxe~ u ^ d(, —y/a<\< 0. (5.15) 

4771 J 

Since dC, = dx + idy, we can split the integral in (15.151) : 

d{ A, a, b ) = <^e~ u{Q dx(Q - ^ <^e~“ (C) dy(Q. (5.16) 

IT IT 

Since the values d( A, a, b), x((), y((), and e~ u ^ are real, the first integral in 
the right hand side of (15.1611 vanishes. So the equality (15.161) takes the form 

d(X,a,b) = (^e~ u(C) dy{(). (5.17) 

IT 

Since the contour Ta is symmetric with respect to the real axis R and the 
function u also is symmetric: u(Q = u((), the equality (15.171) can be reduced 
to the form 

d(X, a, b) = j' e _u(c) dy( C), (5.18) 

rj 

where Tj = Ta n H + is the upper half of the contour Ta- Integrating by parts 
in (15.181) . we obtain 

d(X, a,b) = j* e~ u(a y(() du((), (5.19) 

rj 

(The values y((±(X)) at the end points £+(A), (15. 1311 . of the integration path 
T^ vanish.) 

The differential du(() in (15.191) can be represented as 

MO - ds, (5.20) 

as 

where s is a natural parameteiH on Tj. In other words, the differential du(Q 
can be represented as 

du(() = ^(OdsiQ, (5.21) 

where t(£) is the tangent vector to the curve Tj at the point ( 6 Tj. The 
direction of the vector t(() corresponds to the motion of the point £(s) along 
the path Tj from its left end point £_(A) to the right end point £ + (A). If 
n(() is the vector of the exterior normal to Tj at the point £ e Tj, then 


3 Length of arc. 
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the orientation of the frame (r(£), n(C)) coincides with the orientation of the 
natural frame of R 2 . According the Cauchy-Riemann equations, 


Thus the representation (15.191) can be reduced to the form 


(5.22) 


d( A, a, b) = -L J e"“(« y{ <) J(C) da(C). (5-23) 

rj 


According the assertion 3 of Lemma 15.31 

^(C)>0, VC6T+. (5.24) 

The inequality in (15.241) is strict because ||(£) = I gradu(£)| and the gradient 
gradu(C) of the function v vanishes only at the critical points £+(A) of the 
function v, which are the end points of the integration path Tj. Evidently 
y(() > 0 and e ~ u ^ > 0 at every point £ e Tj. Thus the integrand in (15.231) is 
strictly positive at every point £ e Tj. So the inequality d( A, a,b) > 0 holds. 

□ 


Remark 5.1. The method which we use in the second proof of Theorem [T| is 
the lite version of the method which Herbert Stahl, [Stj, used in his proof of 
the BMV conjecture. 


6. The third proof of Theorem [[] 


For each fixed p, the function cosh (p\rtS + f) is an entire function of the 
variables t, f. Therefore, the Taylor expansion holds 

£) = J] ( 6 - 1 ) 


cosh 


+ i 


0^fc<oo 


where 




df cosh 


+ 


f) 


d£ k 


i«=o 


k = 0,1,2, 


It turns out that for every fixed real p and for every k = 0,1,2, ..., the 
function ipk[t,rj) of the variable t is exponentially convex. We prove this by 
induction in k. Therefore for £ ^ 0, the sum of the series in EH) is an 
exponentially convex function of t. To obtain Theorem [I] we put p = \fa , 
£ = b/a in (16.ID . (For a = 0, the statement of Theorem [T] is trivially true.) 

Our proof of the exponential convexity of the functions 1 [>k(t, if) is based 
on the identity 

= n c<*hA, 

r 11 2 m 

s l^m<oo 

which holds for every (6 C. Substituting the expression 


C = v \A 2 + £ 
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into this identity, we obtain the equality 


sinh (y-\/t 2 + £) 


n c ° sh (^Vi 2 +c)- 

l^m<oo 


Using the equality 

d cosh (j]y/ 1 2 + £) rj sinh (r]y/t 2 + £) 

de = 2 vfr? ’ 

which holds for every t. £, 77 , we obtain the equality 

n («> 

^ l^m<oo 

By successive differentiation the equality (IQ) with respect to £, we 
obtain the equality 

d k+1 cosh (rj^/t 2 + £) _ ? ? 2 ^ f yq d lm cosh + £)\ fR 

^ = UfeUi. ^ > ( } 

where k = 0,1,2,3, ... . In (16.31) . the summation is extended over all se¬ 
quence^] 1 = (h,h, h, • ■ •) of non-negative integers for which |l| = l\ + I 2 + 
I 3 + ... = k. 

The equality (16.311 holds for every t, £, 77 . Restricting this equality to the 
value £ = 0 , we obtain the equality 

i>k+i(t,ri) = yYi( El (6-4) 

" |l| = fc ' l^m<oo ' 

which holds for every t, 77 , and k = 0,1,2,3, ... . In (16.41) . the summation 
is extended over all sequences 1 = ( 1 1 , ^ 2 , ^ 3 , ■ • ■) of non-negative integers for 
which \l\ = h + 1-2 + I 3 + = k. 

Let 77 be an arbitrary real number. By Lemma 1 1.21 the function 


4>o{t, d) = cosh77t (6.5) 

of t is exponentially convex. Moreover, the function tpo (t, is exponen¬ 
tially convex for every m = 1, 2, 3, ... . (The number here plays the same 
role as the number 77 in (16.51) : it is an arbitrary real number.) 

Given k Js 0, assume that all functions ipi(t, ^r) with 0 =% l < k are 
exponentially convex functions of t. Then for each sequence 1 = (R, I 2 ,I 3 , ...) 
with | [| = k, the inequalities 0 < l m ^ k hold. Thus, all the factors ipi m (t, 

which appears in the product J J ipi m I t , ) are exponentially convex 

l^m<00 ' ' 

functions of t. Hence the product itself is an exponentially convex function. 
Finally, the function ipk + i(t,r], (16.411 . which is essentially equal to the sum 


4 For l m = 0, 


d lm cosh (^ y 7 * 2 + ?) 
d£}m 


= f cosh ( jSr \/* 2 + ?) ■ 
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of all such products with |1| = k, is exponentially convex. This finishes the 
proof. □ 

Remark 6.1. Comparing the expansions (I4.2al) and (16.11) . we see that 

<Pk(t,a) = ipk(t,y/a)a~ k , k = 0,1,2,..., let. ( 6 . 6 ) 

As a byproduct of the third proof of Theorem [1] we proved that each of the 
functions (pk(t,a) is exponentially convex. Thus we have given a second proof 
of Theorem [3 


Remark 6.2. Actually we proved more then we formulated in Theorem Q] 
Namely we proved that for any sequence a k (rj) of non-negative numbers the 
sum of series 

s(t) = 2 a k(v)' l Pk(t,ri) (6.7) 

0^fc<oo 


is an exponentially convex function if this series converges for every real t. 
If m is a positive integer and £ ^ 0, then the Taylor expansion 


d m cosh ( g-sjt 2 + £) 

is of the form JO with akiji) 
for k = m, to + 1, to + 2 , ... . 


? Tkhnf^-^ < 6 - 8 > 

m^k<OD v 7 

0 for 0 < k < TO, a k {il) = 7 —rr£ fe_m 

(fc — to) ! 


In particular, for m = 1 the following result holds: 

Theorem 3. For any a > 0 and b > 0, the function ip(t) 
is an exponentially convex function of the variable t. 


sinh (Vat 2 + b ) 
Vat 2 + b 
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